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INTRODUCTION 

Recently the discovery of the charged ressonance state Z(4433) by the Belle Collabora- 
tion [1] has led to the proposal that this state is indeed a radial excitation of a tetratar- 
quark of the form ucdc [2], being the particles X(3872) and X(3876), discovered by 
Belle and Babar, the fundamental tetraquark states ( ucuc and dcdc ). 

Others interpret this state as a scattering ressonance of the D*(2010)Z)(2420) system 
[3, 4] while still others interpret this as a threshold cusp in the same system [5]. 

The Z(4433) was detected by the decay Z(4433) -> tz ± ^, but if the model of the 
tetraquark is correct, the decay Z(4433) — > D ± D° could also be possible. 

SPECTRUM AND WAVEFUNCTIONS 

We start to study the tetraquark system, by using the flux tube model, in wich we have a 
potential given by 

V{r\ , r 2 , r 3 , r 4 ) = V Co uiomb + V Con f 0) 

The term Vcoulomb is the Coulomb part of the interaction corresponding to one gluon 
exchange, and Vconf is the confining part of the potential. 

The Coulomb interaction should be proportional to tftSt, where t\ and ti are the 
Casimir operators of the two particles. Since the two quarks are in an antitriplet state 
we get (t\ + 12) 2 = t\ = ^ = wich gives t\ ■ t% = —2/3 so the quark-quark potential is 

Vm{r) = -\^ (2) 



For the quark-antiquark Coulomb potential, since a quark and two antiquarks form an 
antitriplet, we get a factor —1/3 for the quark-antiquark interaction 

r*M = ~y (3) 

As for the confining part of potential it is given by 

y c onf{r\ , r 2 , r 3 ,r A ) = oL min (r\ , r 2 , r 3 , r 4 ) (4) 

Where L mi „ is the minimum distance that links the four particles as shown in figure 1 
The potential is given by 
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(5) 

We should solve the schroedinger equation for this potential, but because the form of 
the confining part of the potential is cumbersome, we approximate our system by a sys- 
tem of a diquark and a di-antiquark, so that we only have the interaction between quarks 
( antiquarks ) in the diquark ( diantiquark ) and the interaction diquark-diantiquark. Our 
simplified potential is given by 



V (r w ,r c ,rc, tj) = V diq (r u ,r c ) + V diq (r d , r d ) + V diq _ diq (R uc , R c - rf -) (6) 
= a\r u -r c \ + a\rj-r e \ + a\ R uc - R dd \ 

where R uc and R d j are the centers of mass of the uc and dc system. 

, , 2 a s 

V diq (ri , r 2 ) = C diq - - ■ — f + o n - r 2 (7) 

H H 3 |ri — r 2 | 
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V diq - diq (r h r 2 ) =C diq - diq --- ^- + a\ri -r 2 | (8) 

J |Tl — 1*2 1 

So we define the coordinates 

■p _ m u r u +m c r c +m u rj+m c Tc 

K _ 2m u +2m c ^> 

p uc = r u -r c (10) 

Pdc= r d~ r c (ID 

r T = Ruc-Rjc (12) 
Then we could also write the kinetical energy operator as 

f= & vi (13) 

2M 2jUi 2,u 2 2/x v 7 

Writing the wavefunction as 

¥(r M ,r c ,r c -,rj) = ^cM(R)^«c(p MC )V^c(Pjc)^(rr) (14) 



FIGURE 1. String picture for the tetraquark. 



We arrive at three independent Schrodinger equations 

' 2 W + Vy/ = £¥ (15) 

2m 

In this work we only care about s wave particles, so we expand the wavef unctions in the 
s-wave solutions of the harmonic oscilator 

»W=S>.^^«-^ (16) 

n Pi 

We do the same expansion to calculate the meson wavefunctions, wich we use latter 
to calculate the decay width. We use the parameters a = (440MeV) 2 , | a s — 0.27 and 
fix the constant so that we get a mass of 4433MeV for the first radial excitation of 
the tetraquark, and we get the value C = —232%MeV. This is a standard procedure for 
mesons, we choose the potential constant to get the correct masses. The results are given 
in table 1. 

TABLE 1. Masses and constant C 



for different particles 


Particle 


C(MeV) 


M(MeV) 


n 


-1557 


140 




-680 


2980 


J/*p 


-797 


3097 


D 


-976 


1867 


Z(1S) 


-2328 


3913 


Z(25) 


-2328 


4433 


Z(35) 


-2328 


4830 



DECAY WITH MICROSCOPIC QUARK EXCHANGE 

Now we study the decay of the tetraquark system into a two meson system. Note that 
the potential is different in the tetraquark and in the meson- meson channels. We assume 
the "true" potential of the system is the the flip-flop potential 



V F F(r M ,r c ,r d -,rc) =mm(V T ,V„ v ,V DD ) =mm(V T ,V 7t + V v ,V D o + V D +) (17) 



Where Vj is the tetraquark potential, given above, V n , Vy, V D o and Vd+ are the quark- 
antiquark potentials in the respective mesons, wich is of the funnel type 

V i {r i ) = C i -\ C ^ + ar i (18) 

The difference Vff — Vj is the perturbation that causes the decay of the tetraquark into 
two mesons. The full hamiltonean is given by 

H = £7-(p0 +VFF(r„,r c ,rj,r c -) (19) 

We could rewrite the hamiltonean as H = H\ + V\ with V\ = Vff — Vj and V 2 = 
Vff — Vjtiff. Now we calculate the matrix elements 

(T\H\T) =M T + (T\Vi\T) (20) 
Note that we need to redefine the potential so that M T = Mj + (Vi). We have, also 

(tW\h\t) = M^(^|r) + (^|Vi|r) (21) 

(7W\H\tc¥) = M n +My+(TM\V 2 \7i<l!) (22) 
To compute the decay of the tetraquark in two mesons we have to solve the equation 

A*? + <Vi>r M<» e *<0| + v*<x| \/ 1 \ / 1 \ 
M T e\$) + v\x) M K +M V + T+(V 2 ) W J\\x-V) J \\n-Y)J K ) 



where 

1 e*<0| 
e\+) 1 



(24) 



v(r^|^) = ^J rfV^V v ,^«c(p MC )>j c -(pj c -)*0r(r f )V(p ?r ,p v ,,r^)^(p K )Vy{p w ) 

(25) 

e{*n W \§) = g J d 3 p n d 3 p v \ir uc (p uc y\ir d ^Pj d y^T(r t )\ir n (p K )\ir ¥ (p v ) (26) 

The term e\<j>) comes from the nonorthogonality between the tetraquark state and the 
meson-meson state. If the states where orthogonal the S matrix would be the identity. 
We have the equations 

M T -E + (e*M T -e*E)((j)\7t-\i/} + v*(x\n- W) =0 (27) 

(A4-E)e\<j>) + v\x) + (M K + M v + T + {V 2 ) w -E)\7t-xif) = (28) 
Now eliminating \n—\\f) we get 

^ 3 p z*(p)z(p) 



E-M T = -!^- rigmpj 

J (2n) 3 M n +M v + T(v) + (V 2 )K V -E-i£ 



(29) 



where 

X(p) = (*4-E)e(p w \<l>) + v(p KV \x) (30) 

To calculate the decay width, we make the replacement E — > E — ij in the left side of 
the equation and make E ~ in the right side. Also we neglect the term {y 2 )n^ wich 
corresponds to a short range iteraction between the mesons. 

For the terms v(p K y\x) and e(p^|0) we use the eq. 16 for all the wavefunctions, so 
thatv(p^|^) is given by 

HPuc)HpdHPn)HP W )P5{*T)V{p n , Pw ,r T ) 
where the i>, are polynomials. p uc , pj- and r r are given in terms of p K , p v and r nv , by 

P„ c = \(pK~P v )+ r W (32) 
1 

Pdc = 2^Pv~ p ^ +r7tx ^ (33) 

r T = R uc -R dd = p K + p (34) 

m u + m d m u + m c v 

For e{p n xj/\(j>), we just replace the potential by 1. Since the weight of the integral is 
gaussian, we use the Monte Carlo method to evaluate this nine-dimensional integral. So 
we calculate the width for the decay of the Z(4433) ( assuming it is a tetraquark in the 
25 state ) in to ^i//(15) and Ky{2S), and obtain ( preliminary results ) T{ny(\S)) = 
O.lMeV and T(\]/(2S)) = 4.6MeV. 



CONCLUSIONS 

The result is somewhat smaller than the experimental result for the decay of the Z(4433) 
wich has a width of T = 44+}3+} 3 [ MeV [4], but this result is only preliminary. Also, we 
don't include the spin effects, and it is known that the hyperfine splitting effect could 
increase the decay width. 

Even though, we think that our method is quite powerfull and general, and we expect 
to get more accurate and conclusive results using this method. 
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